


o 

(N 



m 
en 

o 

(^ 
-)— > 



o 



^ 



Effect of Depfiasing on Charge-Counting Statistics in Chaotic Cavities 

Ya. M. Blanter''''^, H. Schomerus'', and C. W. J. Beenakker'' 

^Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzerland 

''Laboratoire de Physique Mathematique (CNRS UMR 5825), Universite MontpeUier II, Place 
E. Bataillon, 34095 MontpeUier, France 

"^Instituut-Lorentz, Univcrsitcit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 



Abstract 

We calculate the cumulants of the charge transmitted through a chaotic cavity in the limit that the two openings 
have a large number of scattering channels. The shot noise, which is the second cumulant, is known to be insensi- 
tive to dephasing in this limit. Unexpectedly, the fourth and higher cumulants are found to depend on dephasing: 
A semiclassical theory and a quantum mechanical model with strong dephasing give a different result than a fully 
phase-coherent quantum mechanical theory. 
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1. Introduction 

The theory of shot noise in niesoscopic sys- 
tems has developed along two parallel lines, one 
fully quantum mechanical, the other semiclassical 
(see [y||] for reviews). The semiclassical method 
treats the dynamics of the electrons according to 
classical mechanics, but includes quantum sta- 
tistical effects following from the Pauli exclusion 
principle. The fully quantum mechanical method 
includes interference effects that are ignored in 
the semiclassical approach. Both methods are ex- 
pected to give identical results in the limit that 
the Fermi wavelength goes to zero, or equiva- 
lently, in the limit that the number of scattering 
channels N goes to infinity, corresponding to a 
conductance large compared to the conductance 
quantum e'^/h. 

For example, in a diffusive conductor the shot- 



noise power S 



;(/) in the limit N —> oo 



(with (/) the mean current), and this result has 
been obtained both quantum mechanically (us- 
ing random-matrix theory [pj or Green functions 
El-pl) and semiclassically (using the Boltzmann- 
Langevin equation H). Another example of the 



correspondence principle between quantum me- 
chanics and semiclassics in the large-A'^ limit is 
the shot noise of a double-barrier tunneling diode 

Ml- 

The shot-noise power is the second cumulant 
of the charge transmitted through the conduc- 
tor in a certain time. For the complete counting 
statistics, one needs to know the higher cumulants 
as well. The only existing semiclassical calcula- 
tion of higher cumulants was done for the double- 
barrier tunneling diode |ll| , and was found to 
be in complete agreement with the quantum me- 
chanical theory fl^-flql in the large- A^ limit. 

In this paper we present a semiclassical theory 
for the counting statistics of charge transported 
through a chaotic cavity, using the 'minimal cor- 
relation' approach developed recently by Sukho- 
rukov and one of the authors |lj] . Much to our 
surprise, we do not recover the quantum mechan- 
ical results in the large- A^ limit: the fourth and 
higher cumulants differ. We need to introduce 
dephasing into the quantum mechanical theory 
(by means of a voltage probe p^-|2C|]) to obtain 
agreement with semiclassics. 




Figure 1. Two-terminal chaotic cavity. The left 
and right contacts have cross-section fl^, Qj^, and 
outward normals n^, nj^. 



2. Transmitted charge 

We consider a chaotic cavity (or quantum dot) 
connected to two equilibrium reservoirs by ballis- 
tic contacts (Fig. |]). (Because our results do not 
depend on the dimensionality of the cavity we fo- 
cus on the two-dimensional situation, which sim- 
plifies notation.) A voltage V is applied between 
the reservoirs. The reservoirs are described by 
equilibrium Fermi functions, which at zero tem- 
perature (the case of interest in this paper) have 
the form /l(£;) = 9{eV~E) and fuiE) = 9{-E). 
The cross-sections of the contacts are denoted by 
f2i and Qr. Since we want to exploit the chaotic 
dynamics of the electron motion, we assume that 
there is no direct transmission between the con- 
tacts. In addition, the contacts are narrow as 
compared to the circumference of the cavity, but 
are classical in the sense that the number N^, 
Nji of propagating modes through each contact 
is large. 

We are interested in the statistics of the charge 
transmitted through any of the contacts in a given 
time interval t. 



Qit)= f I{t') 
Jo 



dt' 



(1) 



where /(t) = (/) + 51{t) is the fluctuating cur- 
rent. We take the limit t ^ oo. For definiteness, 
we consider the current through the left contact, 
with the convention that positive current is di- 
rected towards the sample. In the long-time limit 



the charge transmitted through the right contact 
is the same, due to current conservation. 

The average transmitted charge gives the av- 
erage current (/) = \im.t^oot~^{Q{t)), and the 
variance of the transmitted charge gives the zero- 
frequency shot-noise power, 



S = 2 



{SI{Q)SI{t)) dt 



= 2 hm r 



'{{Q')-{Qn 



(2) 



Higher cumulants of Q{t) determine higher-order 
correlation functions of the current. 

Classically, the current is expressed in terms of 
the distribution function /(r, n, E, t), with n the 
direction of momentum. The absolute value of 
the momentum pp niay be regarded as a constant. 
The current is given by 



hit) 



ePF / , 



dE 



/ dn(nL •n)/(x,t). 



(3) 



where ul is the outward normal to the surface of 
the cavity at the location of the left contact (Fig. 
0), yL is a point at the cross-section of the con- 
tact, and we have abbreviated x = {r, n, E}. The 
angular integral is normalized to unity, J dn ~ 1. 
Electrons incident from the left reservoir (n^ • 
n < 0) arc described by the non-fluctuating dis- 
tribution function /i(_E). (We ignore thermal 
fluctuations.) The distribution function of the 
electrons with n^ • n > carries information 
about the chaotic dynamics inside the cavity. On 
average, it only depends on energy, since the coor- 
dinate and angular dependence averages out due 
to the multiple scattering from the surface of the 
cavity. This average distribution function inside 
the cavity is readily found from the condition that 
the current in each energy layer is conserved |17|| , 



(./> = 



NLfL + NnfR 

Nl+Nr ■ 



(4) 



The number of modes N^, Nr is given by Ni^^r = 
PFWL,R/Tf^, with Wl,r the width of the con- 
tacts. The mean transmitted charge M) is thus 



(g) = {iL)t = 



e^Vt NrN, 



R 



h Nl + Ne 



(5) 



This result could also be obtained by series ad- 
dition of the balhstic conductances N^e^ /h and 
N^e^ jh of the two contacts. 

3. Fluctuations of the distribution function 

Fluctuations of the transmitted charge at zero 
temperature are entirely determined by the fluc- 
tuations of the distribution function inside the 
cavity. Following Ref. ||l^, we assume that the 
fluctuations of the distribution function 5/ = 
/ — (/) inside the cavity may be decomposed into 
two parts, 



bj(x,t) = f![x,t)+6!c{E,t). 



(6) 



The function (5/ describes the purely ballistic mo- 
tion and obeys the equation 



(a,+t;fn.V),5/>,t)=0, 



(7) 



with vp the Fermi velocity. This equation is sup- 
plemented by the expression for the equal-time 
correlator, 

(<5/>, t)J/V, t)) = b(x - x'){f) (1 - (/)) , (8) 

where 

6{x - x') = i^-^S{r - r')S{n - n')S{E - E'), (9) 

and v = m/2Trh^ is the density of states. Eq. (||) 
has the same form as in equilibrium, although the 
function (/) describes the non-equilibrium state 
in the cavity. The reason is that the only source 
of noise at equilibrium as well as in the case of a 
deterministic chaotic cavity is the partial occupa- 
tion of states. The uniform fluctuating term Sfc 
accounts for the fact that at long times the mo- 
tion inside the cavity is not purely ballistic, and 
ensures the current conservation at any time in 
any energy interval. 

The fluctuations in the current due to the term 
Sf are given by 



SlLit) 



2 / "^2/1 

inU Jo,, 



dE 



dn{nL-n)Sf{x,t), (10) 



(nfn)>0 

and similarly for SI]i{t). The condition that the 
current at any time in any energy interval is con- 
served can be used to eliminate the fluctuation 



in the current due to the term Sfc- The total 
fluctuation of the current is then given by 



siUt) - 



NRShit) - NL6lR{t) 



Nr+N 



R 



(11) 



Eqs. (Jg) - (g) constitute the minimal correla- 
tion approach, which was shown in Ref. [|l7| to 
agree in the limit N^, Nr 3> 1 with the quantum 
mechanical expression for shot noise in a multi- 
terminal chaotic cavity. We seek to extend this 
approach to higher cumulants of the transmit- 
ted charge. Therefore we have to specify all the 
cumulants of the function Sf at equal time mo- 
ments. Generalizing Eq. m), we write 

{{Sfixi,t)5f{x2,t) ■ ■ ■ Sfixk, t))) = S{X, - X2) 

X S{x, -xs)--- S{x, - xt,){{f\x, t))), (12) 

where ((• • •)) denotes the cumulant. The cumu- 
lant {{f^{x,t))) is calculated later on; it depends 
only on _E, since moments of f{x) are determined 
by (/). Eq. (||) is equivalent to Eq. (|) for fc = 2, 
since ((/2)) ^ (/)(1 - (/)) (see below). 

Taking into account Eq. (|^) and using Eq. ( p^ ) 
as the initial condition, we obtain the expression 
for the cumulant of df at arbitrary times, 

{{6fixi,h)df{x2,t2) ■ ■ ■ dfixk,tk))) 

= i^^-''SiEi-E2)---S{Ek-i-Ek) 
X S{ni - n2) • • • S{nk-i - n^) 
X J [n - r2 - vpuiiti - i2)] • ■ • (5 [rfc_i 
-rk-VFni{tk-i-tk)]{{f'')). (13) 

4. Semiclassical cumulants 

Eqs. d), (0), (0), and ^) may now be used 
to calculate the cumulants of the transmitted 
charge. We note that all cumulants of the type 
{{SIa{ti) ■ ■ ■ (5/,,(ifc))) vanish due to the combina- 
tion of delta functions in Eq. ( [l3|) , unless all the 
subscripts a, . . . , 7y are equal to L or all are equal 
to R. The two non-zero cumulants are 



{{sL{h)sL{t2)---sUtk))) = {-i 



.k^'N^ 



X / dE{{f''{E)))5{h - t2) ■ --Sitk-i - h), 



(14) 



a:^L,R, fc > 2. 



Consequently, for the cumulant of the transmit- 
ted charge we obtain 



{{Q')) = 



~~h {Nl + NrY 

X [dE{{f\E))). 



(15) 



To complete the calculation, we must compute 
the cumulant ((/'')). The fluctuating distribution 
function assumes only the values and 1, thus 
if^) = (/) for fc > 1. The characteristic function 



X{p) = (exp(p/)> - 1 + (/) [exp(p) - 1] 



(16) 



generates the cumulants as coefficients in a series 
expansion. 



i-x(p) = Et7((/')) 



k=l 



(17) 



Substituting the average distribution function 
(y), integrating over energy, and taking into ac- 
count Eq. (H), we obtain the expression for the 
cumulants of the transmitted charge (fc > 2), 



{{Q'))^e^~'{Q)S,^'^''^^ ^^'^«"' 



fe-i 



{Nl + Nr) 



where the coefficients Si- are defined as 



\n[NR + NLeMp)] =J2Ti^'' 



fe=0 



(18) 



(19) 



Expression (n8[) may be simplified in the sym- 
metric case Nl = Nr, when 



_ 22' - 1 
S21 — — ^1 — B21, 



S21+1 = 0, ; > 1, 



(20) 



where Bn are Bernoulli numbers. In this case we 
have (for / > 1) 

22' -1 






((Q"» = '-'''(Q) 1221^1 



B 



21, 



((Q-+^» 



0. 



(21) 



In particular, {(Q"^)) = e(Q)/4, which is the 1/4- 
shot noise suppression in a symmetric chaotic cav- 
ity [ pl| . The next non- vanishing cumulant is neg- 
ative, ((Q4))=-e3(Q)/32. 



For non-equivalent contacts Eq. { {wj) yields 



{{Q'))=e{Q) 



NlNr 



(Nl + Nr^ 



in agreement with Refs. 
mulants are 



(22) 
The next two cu- 



l4\\ _ „3 



NLNRiNL - Nr)^ 



[Nl + Nr)^ ' 
((Q"» = e-\Q)NLNR 

^^^ {Nl-mLNR + Nl){Nl + Nl,) 
{Nl + NrY 



(23) 



(24) 



5. Phase-coherent quantum mechanical 
cumulants 

We produced the result ( |lq ) by a semiclassical 
method, generalizing the minimal correlation ap- 
proach to higher cumulants. Let us compare it 
to the phase-coherent quantum mechanical the- 
ory in the large- TV limit. 

We use the relationship derived by Lee, Levi- 
tov, and Yakovets [|l5[ between the characteristic 
function of the transmitted charge and the trans- 
mission eigenvalues Tj, 

{expiQp/e)) = n [1 + Tj (exp(p) - 1)]^^*/'' . 

(25) 

This is the expression for a particular cavity, so 
we still have to perform an ensemble average. 
The transmission eigenvalues have density p{T) 
in the ensemble of chaotic cavities. The ensemble- 
averaged cumulants of the transmitted charge fol- 
low from 

k=l 



eVt 
IT 



f dTp(T) In [1 + T(exp(p)-1)], 
Jo 



(26) 



The density of transmission eigenvalues for Nl, 
Nr > 1 has the form flpl 



P{T) 



VWNr 

ttT 



T 



1 -T 



[Nl - Nr) ' 
4.NlNr 



2\ 1/2 



(27) 




Figure 2. Chaotic cavity with a fictitious lead 
(supporting N^ propagating modes) that models 
strong uniform dephasing. 



It vanishes for T < [1 + ANlNr{Nl - Nr)-"^]-^. 
In the symmetric case N]^ — Nr — N one has 
simply MM 



P{T) = 



N 



1 



T)' 



(28) 



4(Q)ln[i + iexp(p/2)],(29) 



vr v/T(l- 
and hence 

oo ^ 
k=l 

which is equivalent to Eqs. (nSl) and (E3). Thus, 
for the symmetric case the results for the charge 
counting statistics, obtained semiclassically and 
quantum mechanically, are identical. 

This correspondence does not extend to the 
more general case Nl ^ Nr. From Eqs. (|2q ) and 
(27) we reproduce the semiclassical expressions 
@ and dl) for {{Q^)) and ((Q^)). However, 
the fourth cumulant, 



((Q"*)) = e''{Q)NLNR{NL + Nr) 



x{N 



^ - snInr 



UNJNl 



SiViiVl + TV*), 
(30) 



is difi^erent from Eq. (P4| 



6. EfTect of dephasing on the quantum me- 
chanical cumulants 



Apparently, the difference between Eqs. (24) 
and (BG) is due to the lack of phase coherence in 



the semiclassical minimal correlation approach. 
To strengthen this explanation, we incorporate 
into the quantum mechanical theory a dephas- 
ing mechanism via an additional fictitious lead 
attached to the cavity (see Fig. ||). Following 
Refs. ||l^ and ^^, we assume that this lead has 
a (fiuctuating) distribution function which is de- 
termined from the condition that no current fiows 
into the reservoir at every energy and every in- 
stant of time. An electron absorbed by this lead 
is immediately re-injected back at the same en- 
ergy, without any memory of the phase. Thus, 
this fictitious lead introduces dephasing but not 
inelastic scattering. For homogeneous and com- 
plete dephasing, we take the number of scatter- 
ing channels iV^ in the dephasing lead very large, 

N^:$>Nl, Nr. 

Our starting point is the expression for the cur- 
rent operator in the lead a {a — L, R,(j)) in terms 
of creation/annihilation operators |^,||, 

X al„,{E)A:;r{o',E,E')asn{E'). (31) 

The sum over 7 and S is over the leads L, R, and 
0, the sum over m, n is over the mode indices. 
The current matrix A is expressed in terms of 
the scattering matrices s. 



A™,"(a,i?,£;') 



SmnSajSaS - si'^"{E)sZ\E'), 



(32) 

and the creation and annihilation operators obey 
the fermion anticommutation relation 

al^{E)apn{E') + ap,,(E')al^{E) 

= 5c.p5,nn5{E-E'). (33) 

The expectation value 

{al„,{E)apr^{E')) = 5o.p5.^n6{E - E')U{E) (34) 

is given by the average distribution function 
Ja{E) in reservoir a. From the condition that 
the average (energy resolved) current {I^{E)) 
through the dephasing lead vanishes at every en- 
ergy, one finds the average distribution function 



UiE) 



NlIl + NrJr 

Nl + Nr ' 



(35) 



which is identical to the average semiclassical dis- 
tribution function (jj) inside the cavity. 

The distribution function of the dephasing lead 
fluctuates in time. These fluctuations, which can 
also be found from the condition that the current 
through this lead vanishes, modify the fluctua- 
tions of the current at the leads L and R |20l , 



A/L(t) = <5/L(i)- 

A4(i) = Sidt) 



Nl 



Nl + Nr 
Nr 



Nl + Nr 



suit), 



(36) 



where the intrinsic fluctuations 5Ia{t) = Ia{t) — 
(Lit)) are described by Eq. @. 

We now calculate the fourth cumulant of the 
charge transmitted through the left lead, 



((Q'» = (((£rfi'A/L(i')) »• 



(37) 



Explicitly, we have 

((Q'»=s 



■'LLLL 



ml 



^Nl „ 

^LLL4> 



(Nl + Nr) 

Nt 



Nl + Nr 



ANf 



[Nl + Nr)^ 



{Nl + Nr)^ 
where we have defined the cumulant 



(38) 



(39) 



in terms of the current operators (|3l|). It re- 
mains to calculate the cumulant of eight cre- 
ation/annihilation operators, by taking into ac- 
count all possible pairings between the operators 
at^-rui ^^'^ ^SiJii {i = 1, 2, 3, 4) which couple all 
the four current operators together. For N^p ^ 
Nl , Nr ^ 1 the leading contribution to the cu- 
mulant comes from the terms with "fi = Si = (j). 
Summing all possible pairings, we obtain 



eH 



-"010^20^30^4 



Tr [A^^{ai) ■ ■ ■ A^^{ai)\ / dE 



-X U{1~ U) 



id, 



.{i-f^) + [i-f^y 



(40) 



where the trace is taken over the mode indices. 
In the same leading order, we neglect all traces of 
the type Tr [s^^ • • • s^^] (a ... 77 = L,R), unless 
all indices a . . .rj are the same. In this case 



Tr 



Lrj)^Lit> 



^Nr 



Tr 



^Nn 



(41) 



Summing all contributions in Eq. (pTJ), we recover 
the result ( |24| ) of the semiclassical theory. 

7. Conclusions 

To summarize, we have studied the charge 
counting statistics in chaotic cavities using three 
different approaches: (i) a fully coherent quantum 
mechanical theory; (ii) a dephasing-lead model 
that is also quantum mechanical but phenomeno- 
logically introduces uniform and strong dephasing 
inside the cavity; and (iii) the semiclassical min- 
imal correlation approach. All three approaches 
give the same results if the two openings in the 
cavity have the same (large) number of scattering 
channels. The cumulants for this symmetric case 
are given by Eq. (|2l|). In the asymmetric case of 
two different openings the first three cumulants 
are also the same in the three approaches, but 
the fourth cumulant is not: For approach (i) it is 
given by Eq. (^0|), for approaches (ii) and (iii) by 
Eq. (H). 

We conclude that the counting statistics of a 
chaotic cavity is sensitive to dephasing even in 
the limit of a large number of scattering chan- 
nels. This is an unexpected conclusion since the 
shot noise is not affected by dephasing. What is 
lacking is an understanding in physical terms for 
why the high-order cumulants of the transmitted 
charge respond differently to dephasing than the 
low-order cumulants. 
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